Salt-induced counterion-mobility anomaly in polyelectrolyte electrophoresis 
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We study the electrokinetics of a single polyelectrolyte chain in salt solution using hydrodynamic 
simulations. The salt-dependent chain mobility compares well with experimental DNA data. The 
mobility of condensed counterions exhibits a salt-dependent change of sign, an anomaly that is also 
reflected in the counterion excess conductivity. Using Green's function techniques this anomaly is 
explained by electrostatic screening of the hydrodynamic interactions between chain and counterions. 



Polyelectrolytes (PEs) are macromolecules with ioniz- 
able groups that dissociate in aqueous solution and thus 
give rise to a charged PE backbone and a diffusely bound 
cloud of neutralizing counterions Numerous applica- 
tions in chemical, biological, and medical engineering rely 
on the response of PEs to externally applied electric fields 
(E- fields), determined by a balance of electrostatic and 
hydrodynamic effects and controlled by various factors 
such as salt concentration, PE charge density, etc. 
The simplest scenario providing a basic testing ground 
for our understanding of PE dynamics in the dilute limit 
is free-solution electrophoresis, where a single PE chain 
is subject to a homogeneous static E- field 0, 0, 

Previous theoretical approaches combined mean-field 
electrostatics with low Reynolds number hydrodynam- 
ics. Solutions of the electrokinetic equations were ob- 
tained numerically [91 or analytically using counterion- 
condensation theory [f] and account for the experimen- 
tally measured salt dependent electrophoretic mobilities 
of biopolymers such as DNA or synthetic PEs. Counte- 
rions in the immediate vicinity of the PE chain were as- 
sumed to stick to and move along with the PE under the 
action of the applied E-field. This assumption becomes 
crucial for the conductivity of PE solutions, and indeed 
inconsistencies between experimental mobility and con- 
ductivity studies are documented in literature, pointing 
to some basic riddles in the coupling of PE and counte- 
rion dynamics in E-fields 8]. Pioneering explicit- water 
all-atomistic simulations of PEs in E-fields have been 
performed [9|]. Due to the immense computational de- 
mand they are restricted to elevated field strengths, short 
PEs, and short simulation times. Implicit-solvent simula- 
tions have quite recently addressed the molecular- weight- 
dependent PE mobility in the salt-free case 

[13, El and 

yielded good agreement with experiments. 

In the present paper we use coarse-grained implicit- 
solvent hydrodynamic simulations [l^ and study the 
salt-dependent electrophoretic response of a single PE. 
By replicating the PE periodically we eliminate finite- 
chain-length effects. We concentrate on the salt- 
dependent interplay of PE versus counterion mobility 
in the infinite chain limit and show that the condensed 
counterion mobility changes sign as a function of salt con- 



centration. For low salt, counterions stick to the PE and 
move along in the E-field in agreement with the canonic 
viewpoint. For high salt, on the other hand, the motion 
decouples and counterions move opposite to the PE. This 
anomaly is captured by an analytic theory developed here 
for weakly charged chains based on the electrostatically 
screened hydrodynamic interaction tensor. For DNA our 
simulations reproduce experimental salt-dependent mo- 
bilities without fitting parameters and predict an exper- 
imentally measurable anomaly of the counterion excess 
conductivity. The counterion anomaly is also directly 
accessible by NMR experiments [13] or PE conductivity 
studies in nanopores or nanochannels [l3 |. 

In our hydrodynamic simulations we consider a PE 
consisting of charged beads together with neutralizing 
counterions and added symmetric salt, cf. Fig. [T] The 
vertical box height H and lateral width D are fluctu- 
ating while keeping the volume HD^ and thus the con- 
centration of monomers Cm, neutralizing counterions Cnct 
and salt ion pairs Cg fixed. Periodic boundary conditions 
along the vertical axis are implemented by coupling the 
box height H to the vertical PE extension. All particle 
positions evolve according to the position Langevin 
equation, r,{t) = - J2j My • VrjU{t) + ^i{t). The ther- 
mal coupling is modeled by a Gaussian white noise 
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FIG. 1: Simulation cell for a DNA segment with counterions 
(dark grey) and colons (light grey). Periodic boundary con- 
ditions are applied along the projected end-to-end distance 
H of the DNA segment. The external electric field is applied 
either parallel (-B||) or perpendicular {E±) to the PE axis. 
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FIG. 2: Hydrodynamic simulations of DNA in aqueous NaCl solution of various ionic strengths / at 20 °C. a) Electrophoretic 
DNA mobility fj,ui as obtained from simulations (filled symbols) and experiments (open symbols) as a function of /, compared to 
theories by Stigter (dashed line) and Manning Q (dotted line), b) Mobility of neutralizing counterions finct, and condensed 
counterions Hcct- c) Counterion excess conductivity And according to Eq. ([T]). 



with (|,(t)) = and {^dt)ij{t')) = 2kBTM,^S{t-t') 
according to the fluctuation-dissipation theorem. Hy- 
drodynamic interactions are included via the Rotne- 
Prager-Yamakawa mobility tensor M^- which ac- 

counts for finite hydrodynamic particle radii (a^ = 
Q-m , Q-ct , Q-co for monomers, counterions, and colons). The 
interaction potential U — ?7l.i + Uq + Us + Ucxt con- 
sists of: i) A truncated, shifted Lennard- Jones poten- 
tial, f/Lj/kBT = eEfe) [KAy)''-2K/r,,)6 + l] 
for Tij < (Jij between ions and monomers that pre- 
vents electrostatic collapse of opposite charges, where 
Tij = {fi — fjl is the distance between particles i and 
j and aij — (oj -I- aj)/2 and e define the soft-core dis- 
tance and repulsion strength, ii) An unscreened Coulomb 
potential Uq/^bT = ^B^^iij) lilj /''^iji where qi denotes 
particle valency (g^ — q^, q^, (?co for monomers, counteri- 
ons, and colons) and 1-q = ^ j A.'Kti-tQ&^T is the Bjerrum 
distance at which two unit charges interact with ther- 
mal energy keT (Ib = 7.21 A in water at 20 °C). iii) 
A harmonic potential, t/s/keT — {K /2)Y^i^^^^{r.ij — b)^, 
which acts between adjacent monomers only and ensures 
chain connectivity, iv) The external electric potential, 
Ucxt/kBT — — X)i(9»e/kBT) E-r^, with the electric field 
directed either parallel (£^||) or perpendicular {E±) to 
the PE axis. Periodic boundary conditions along the PE 
axis are implemented by a one-dimensional resummation 
of the Coulomb interactions [15,] ; the lateral and all hy- 
drodynamic interactions are treated using the minimum 
image convention. Consequential finite-size effects are 
discussed in the supplementary information [16]. The 
PE electrophoretic mobility /im = (i;,n) /E follows from 
the average monomer velocity along the E-field direction. 
In the absence of curvature, inter-chain and end effects 
(i.e. for high enough salt concentrations) and if orien- 
tation effects are negligible {i.e. for small E-fields), /im 
follows from the parallel and perpendicular mobilities as 
Mm = (mIIi + 2/ijJj)/3. In the simulations we accordingly 
determine /^m and separately by applying E-fields par- 



allel and perpendicular to the PE axis and measuring the 
corresponding velocities. Possible non-linear effects have 
been carefully checked [l^]. The ionic strength includes 
contributions from the neutralizing counterions and is 
defined as / = (cnct^ct + Cs^ct + Cs9co)/2- 

In order to model DNA in aqueous NaCl solution at 
20 °C we use Stokes radii of Na+ and CI" as act = 1.84 A 
and Onn = 1.29 A as obtained from limiting conductiv- 
ities [l3|, an estimate of = 10.47 A for the DNA 
radius and valencies get — 1, 9co = —1 and g„i — —6. 
The choice of monomer separation b — Om ensures a 
linear charge density of qm/b ~ 0.57A . Although 
no bending rigidity is present in the model, the seg- 
ment is sufficiently straight due to electrostatic repul- 
sions, as appropriate for DNA (cf. Fig. [ij. The simu- 
lation cell comprises 10 DNA monomers, 60 neutraliz- 
ing counterions and 24 salt pairs. The ionic strength 
is varied over the range / = 19 — 468 mM by adjust- 
ing the cell width D. The field strengths applied are 
^li = 27.5 X 10^ V/m and Ej_ = 5.5 x 10*^ V/m. We use 
e = 5 for the LJ strength, K/a^^ — 100 for the bond stiff- 
ness, and 1] = 1.003 x lO^'^ Pas for the viscosity of water. 
The Langevin time-step is 0.06 — 0.12 ps and simulations 
are typically run for 0.3 — 4.1 |j,s. 

In Fig. [2^ we plot the DNA electrophoretic mobility 
Mm as a function of the ionic strength, /, together with 
experimental data for long DNA from Refs. [3, |3]. Not- 
ing that there arc no free fitting parameters and given 
the substantial scatter in the experimental data, we con- 
clude that our coarse-grained DNA model is quite accu- 
rate. The mobility /Xm decreases with increasing / which 
will be rationalized in terms of hydrodynamic screening 
effects below. We additionally show theoretical results 
from Stigter and Manning [3]. 

Theoretically, only little attention has been paid to E- 
field-induced counterion dynamics in PE solutions. In 
this context the phenomenon of counterion condensation 
at highly charged PEs that are characterized by a Man- 
ning parameter = l9ct9m^B/fel > 1 has to be taken 
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into account. For highly charged PEs such as DNA 
(^M = 4.17) electrostatic attraction of counterions to- 
wards the PE overcomes entropic repulsion giving rise to 
increased accumulation of counterions in the very vicin- 
ity of the PE [3, [ill- In particular, the assumption 
that condensed counterions stick to the PE 0, 0] has not 
been scrutinized, despite experimental evidence that con- 
densed counterions are not immobilized on the PE sur- 
face [3]. In Fig. [2b we show the electrophoretic mobihty 
of two counterion ensembles, first condensed counterions 
within a distance = Om -H4act = 17.8 A from the DNA 
axis (/Xcct) and secondly the set of counterions closest to 
the DNA axis that neutralize the DNA charge (/inct)- 
At low ionic strength the hydrodynamic drag exerted by 
the DNA on the counterions exceeds the external elec- 
tric force and the mobility for both sets is negative, i.e. 
the counterions are dragged along by the PE. At high 
ionic strength the hydrodynamic interactions arc suffi- 
ciently screened so that the electric field dominates and 
the counterions move opposite to the DNA. In fact, a 
salt and PE charge density dependent sign reversal of 
the electrophoretic counterion mobility has been inferred 
from transference experiments some time ago jl9| . Di- 
rect measurements of counterion electrophoretic mobili- 
ties can in principle be performed with pulsed field gra- 
dient NMR pj]. 

We define the excess contribution of the counterions to 
the conductivity of a PE solution as 

^nct = (o- - 9meC„i/i,n - Cr°) / Cnct (1) 

where a and denote the specific conductivities of 
the salt solution with and without the PE chain, re- 
spectively. In our simulations, cr results from the sepa- 
rate electrophoretic contributions as a/e — qctCnctMnct — 
gmCmMm + gctCs^ct " gcoCs^co, whilc the purc electrolyte 
conductivity ctq is obtained from separate simulations as 
a° = (jctCsMct - QcoCsfJ-co- As seen in Fig. [2fc, the coun- 
terion excess conductivity A^ct increases with increasing 
salt concentration and changes sign, and thus directly 
reflects the salt-dependent counterion mobility anomaly 
for the experimentally easily accessible conductivity. 

To gain further insight, we now shift to weakly charged 
PEs (Manning parameter < 1), where the ion dis- 
tribution around a PE is correctly described by linear 
Debye-Hiickel (DH) theory and the electrophoretic mo- 
bilities of PE and ions can be constructed using Green's 
functions Q . The DH ionic charge distribution around a 
sphere of radius a and surface charge density qe/Aira'^ is 
for Kfl < 1 given by n{r) — —qeK'^e~'^^^~'^^ /['i'Kr{l + na)], 
where — (Stt^b^)^^^^ is the DH screening length. 
On the Stokes level, the solvent flow field u induced by 
an external electric field E acting on the ionic charge dis- 
tribution n{r) is u{r) — G{r) qeE, where the screened 



hydrodynamic Green's function Gap{r) reads [l3| 

G^p{r) = aI^S^p - + 25^, (2) 

where a,/? — x^y,z and i? = cT'^''^'^'^^ /[ATTrir{l + Ka)\ and 
A = B[l + nr + K^r"^ - e'-'^''-"-\l + Ka + kV/3)]/(kV). 
In the limit of zero salt k ^ 0, the Stokes solution for a 
translating sphere is recovered [13] ■ For vanishing radius 
a 0, Eq. ^ reduces to a previously derived expres- 
sion [21|. Noting that Eq. ([2]) fulfills the no-slip condi- 
tion on the sphere's surface, its electrophoretic mobility 
follows from /ig = qe Gxx [f = a) as 

J^ = ^, (3) 
ge/io 1 + KG 

which is the classical result derived by Debye and 
Hiickel . Here ^J.a = I /GTrrya is the Stokes mobility. 

To leading order in a/r the electrophoretic coupling 
matrix /i"/(T") {a, (3 — x,y,z) between two charged 
spheres is obtained via a multipole expansion [20| as 

f,ffir) = qe(l + ^Vl^G^pir). (4) 

Approximating the PE as a straight chain of charged 
spheres at spacing b oriented along the z-axis, the 
PE mobilities follow by superposition as /i^ = Ms + 

2 Ejli ^^ss{r± = 0, bj) and fiL ^ ^ls + 2 J^'^^i ^^Isi^^ = 
0, bj). Here r± denotes the lateral distance from the PE 
axis. For the orientationally averaged mobility /i,n = 
{pm + 2/i^)/3, we obtain the closed-form expression 

^ = ^ - ^^±^ In (1 - e-^") . (5) 
ge/io ?e/io 36(1 -I- na) 

In the limit of low screening, kq 0, Eq. (O decays log- 
arithmically with increasing ionic strength as /^m/^e/io = 
—2{a/b) ln(K&), in accord with previous results for weakly 
charged PEs |l|, 0| • In the same fashion the perpendic- 
ular and parallel distance-dependent ion mobilities fol- 
low as A*™/ct(''-L, z) = fis± Ejtr^oo -"Sf (''-L, z + bj) and 

f^L/ct (''-L , 2) = Ais ± J2t=-oo l^ls ir±,z + bj), respectively; 
the plus/minus sign applies to coions/counterions. 

In Fig. [3] we compare the foregoing theoretically 
predicted electrophoretic mobilities of monomers and 
ions (obtained by summing over contributions from 23 
spheres) to the hydrodynamic simulations of a weakly 
charged PE with Manning parameter = 0.4 in a field 
of strength aqeE /\-bT = 0.2. The simulation cell com- 
prises 24 PE monomers, 24 neutralizing counterions and 
24 salt pairs with equal radii a,n = Oco = Oct = a, 
valencies q^i = qco — —qct = q and monomer spac- 
ing b ~ 2a. For intrinsically flexible PEs, the straight 
PE conformation in our simulations and theory is real- 
istic only for low enough salt concentration as long as 
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FIG. 3: Results for a weakly charged PE with Manning parameter ~ |gctqm^B/&| ~ 0.4. a) Colon mobility ^ico{r±,0) 
for various fixed distances r± from the PE as a function of (sa)^. The hydrodynamlc drag exerted by the PE on the colons 
Increases their mobility as compared to the case without PE, Eq. (|3} {r_L/a oo, solid line). The latter compares well with 
hydrodynamlc simulations of a simple salt solution (filled symbols) . b) Counterlon mobility /ict {r± , 0) at distance r± from 
the PE which exhibits a sign change. For small r± and low salt, counterlons are dragged along with the PE (/ict(r_L,0) < 0). 
c) Comparison of theoretical predictions (solid lines) and hydrodynamlc simulations (filled symbols) for the PE mobility ^m, 
neutralizing counterlon mobility /inct and condensed counterlon mobility /icct . The condensed counterlon mobility changes sign. 
We also show the parallel and perpendicular PE mobilities /xjL (dotted line) and /x^ (short dashed line). 



the effective persistence length is larger than the screen- 
ing length. In Fig. [3^,b we show the orientationally av- 
eraged colon and counterlon mobilities /^co/ct('''±7 0) = 

[2Mm/ct(^-L' 0) + A*!!o/ct(''-L' for fixed vertical coordi- 
nate z = and various fixed distances r± from the PE as 
a function of the rescaled salt concentration (ko)^ Cs- 
The mobilities of colons are Increased and those of coun- 
terlons are decreased by the presence of the PE. This 
entraining effect is larger for smaller salt concentration 
and smaller r_L. The ion mobilities for r± = oo reflect 
pure electrolyte friction effects and in Fig. [3^ compare 
very well with simulation results for a simple salt solu- 
tion. In Fig. [3fc we compare analytical predictions for 
the PE mobility the neutralizing counterlon mobil- 
ity fJ-nct, and the condensed counterlon mobility /icct (ob- 
tained from counterlons within a shell of = 5a around 
the PE) with the simulations. Here finct and fleet are 
obtained from /ict (f^ > 0) by spatially averaging over the 
DH counterlon distribution around a straight chain of 
charged spheres at fixed vertical coordinate z — 0. With 
increasing salt concentration, ^cct changes its sign, sim- 
ilar to the DNA results (cf. Fig. This shows that 
the salt-induced counterion mobility anomaly is not re- 
stricted to the non-linear regime and is fully explained by 
screening effects of the hydrodynamlc coupling tensor. 

Our simulation method neglects local solvation and 
DNA structural effects. The good agreement between ex- 
perimental and simulation results could Imply that those 
effects are of minor Importance for the electroklnetlc be- 
havior. Nevertheless, an extension of the model to more 
realistic charge distributions is planned. Likewise, the 
analytic Green's function approach will be generalized to 
Include relaxation effects in addition to retardation. 

This work was supported by the German Excellence 
Initiative via the Nanosystems Initiative Munich (NIM) . 
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